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ABSTRACT 

The triple nucleus of M31 consists of a population of old red stars in an eccentric disk (PI and P2) 
and another population of younger A stars in a circular disk (P3) around M31's central supermassivc 
black hole (SMBH). We argue that PI and P2 determine the maximal radial extent of the younger A 
star population and provide the gas that fueled the starburst that generated P3. The eccentric stellar 
disk creates an m — 1 non-axisymmetric perturbation to the potential. This perturbed potential 
drives gas into the inner parsec around the SMBH, if the pattern speed of the eccentric stellar disk 
is fl p < 3 — lOkms -1 pc -1 . We show that stellar mass loss from PI and P2 is sufficient to create a 
gravitationally unstable gaseous disk of ~ 1O 5 M every 0.1 — 1 Gyrs, consistent with the 200 Myr age 
of P3. Similar processes may act in other systems to produce very compact nuclear starbursts. 
Subject headings: galaxies: individual (M31) - galaxies: nuclei - galaxies: starburst 
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1. INTRODUCTION 

In balloon-borne experiments, Light et al. (1974) dis- 
covered that the nucleus of M31 is asymmetric. Obser- 
vations with the Hubble Space Telescope (HST) resolved 
the nucleus into two components (Lauer et al. 1993). 
The two nuclei, denoted PI and P2, have an angular 
separation of 0".5 («s 2 pc at the distance of M31; Ben- 
der et al. 2005, hereafter B05). P2 is located near the 
dynamical center while PI, the brighter nucleus, is offset. 
PI and P2 have a total luminosity of ks 3 x 10 6 L Q . For 
a mass-to-light ratio of 5.7 appropriate for a bulge pop- 
ulation (Tremaine 1995, hereafter T95), the total stellar 
mass is rs 2 x 10 7 M Q . 

PI and P2 are unlikely to be separate stellar clusters 
because the time for dynamical friction to merge such 
clusters would be short (~ 10 6 yrs if PI is a ~ 10 7 M© 
cluster; T95). Instead, PI and P2 are best modeled as 
a single eccentric stellar disk, as originally proposed by 
T95 (see also Peiris & Tremaine 2003, hereafter PT03). 
Disk stars slow near apocenter, giving rise to PI. Stars 
near pericenter, and at small disk radii, give rise to P2. 
The central supermassive black hole (SMBH) sits within 
P2. 

Spectroscopically, PI and P2 are nearly identical, con- 
sistent with the eccentric stellar disk model. However, P2 
is bluer. Neito et al. (1986) showed that P2 is brighter 
than PI at 3750 A. King et al. (1995) showed that P2 is 
brighter than PI in the ultraviolet. Though this differ- 
ence was initially attributed to an active galactic nucleus 
(AGN), recent HST spectroscopy of the nuclear region 
of M31 has uncovered a younger population of 200-Myr 
old A stars, embedded in P2 (B05). This population, 
named P3, appears to be a disk of stars with a mass of 
~ 4200M© and a maximal radial extent of « 1 pc (B05; 
Bender priv. communication) that surrounds the central 
SMBH and lies in the same plane as the P1/P2 disk. Us- 
ing the line-of-sight velocities (« 1000 km s -1 ) measured 

1 Astronomy Department and Theoretical Astrophysics Cen- 
ter, 601 Campbell Hall, University of California, Berkeley, CA 
94720; pchang@astro.berkeley.edu, rmurray@astro.berkeley.edu, 
echiang@astro.berkeley.edu, eliot@astro.berkeley.edu 

2 Miller Institute for Basic Research 



for P3, B05 estimate that the mass of the SMBH in M31 
is 1.1 - 2.1 x 10 s M Q . 

P3 is a stellar population distinct from PI and P2 
(B05). P3 is composed of A stars while PI and P2 are 
typical old red bulge stars. The relative youth of these 
A stars and their proximity to the central SMBH make 
P3 analogous to the young stars in our Galactic Center 
(GC). Like the young stars in the GC, P3 must either 
form in situ or migrate in from larger radii. Migration is 
less likely in M31 than in the GC as the progenitor cluster 
would be disrupted at greater distances from the more 
massive SMBH in M31. In situ formation more naturally 
explains the masses of these central star-forming regions 
through Toomrc stability arguments (see § 3). However, 
it is less clear what sets the radial extents (r < 1 pc for 
M31; r < 0.4 pc for the GC) and ages (« 200 Myr for 
M31; « 10 Myr for the GC) of these nuclear starbursts. 

In this paper, we address these questions by demon- 
strating that the eccentric stellar disk of M31 fixes both 
the radial extent and the timcscale for the starburst 
that generated P3. In § [2 we argue that the non- 
axisymmetric potential from the eccentric stellar disk 
limits non-intersecting gas orbits to a limited family 
around the central SMBH. The physics is similiar to what 
sets the maximum sizes of accretion disks in Roche-lobe 
filling binaries (Paczynski 1977; Papaloizou & Pringle 
1977). We present numerical and analytic calculations 
describing how non-intersecting gas orbits are only al- 
lowed for r < 1 pc if the pattern speed of the P1/P2 disk 
is < 3 — 10 km s _1 pc -1 . This naturally explains the size 
of P3. We then argue in § [3] that stellar mass loss from 
the P1/P2 disk is sufficient to supply the gas needed to 
form P3. We estimate the mass required to trigger a 
starburst and the timescale to build up that mass, and 
show that these are consistent with the mass and age 
of P3. Finally, we conclude in § 01 drawing attention 
to a few predictions of our model and arguing that this 
mechanism may be common in galactic nuclei. 

2. CLOSED GAS ORBITS IN AN ECCENTRIC STELLAR 

DISK 

In the limit that gas has zero pressure, gas follows test 
particle orbits that are simply closed and non-crossing 
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(Paczynski 1977). Gas not in these orbits will shock 
against gas in neighboring orbits (i.e., crossing) or itself 
(i.e., not simply closed). These shocks dissipate orbital 
energy, driving gas to lower energy orbits. Paczynski 
(1977) applied this principle to solve for the maximum 
size of a gaseous accretion disk in a Roche-filling binary 
star system. Test particle orbits that are close to the 
accretor, where the potential is nearly Keplerian, can be 
nearly circular and non-crossing. Farther from the ac- 
cretor, the non-axisymmetric component of the potential 
(due to the donor star) becomes larger until test particle 
orbits are forced to cross either their neighbors or them- 
selves. Therefore there exists a maximum radius for gas 
orbits in the vicinity of the accretor: the tidal truncation 
radius, i? t (Papaloizou & Pringle 1977). Gas outside Rt 
will be driven toward it through dissipative processes, 
while gas inside Rt will occupy the allowed orbits, form- 
ing a gaseous accretion disk. Paczynski (1977) showed 
that only one family of orbits is possible for any given 
binary system with a specified mass ratio. These results 
were later confirmed by numerical simulations of close 
binary systems by Ichikawa & Osaki (1994). 

By analogy, PI and P2 add a non-axisymmetric com- 
ponent to the point mass potential of the central SMBH 
in M31. Thus, there should also be an R t inside of which 
a gaseous accretion disk can exist around the SMBH. 
While this situation is similar to that of a close binary 
star system, there are two differences. First, the per- 
turbation to the potential in M31 is given by the ec- 
centric stellar disk. Second, whereas the pattern speed 
of the perturbation potential in a binary star system is 
prescribed by the binary orbital frequency, the pattern 
speed of the P1/P2 disk (i.e., its apsidal precession fre- 
quency) is uncertain (see Appendix A for estimates in 
the literature). 

Since the gas mass required to form P3 (A/ gas < 
10 5 M©; see § [3]) is much smaller than the mass of the 
P1/P2 disk and the central SMBH, and since we assume 
that Toomre's Q > 1 for the present discussion (§[3]), we 
neglect gas self-gravity. Moreover, since the characteris- 
tic temperature of the gas is T ~ 30 K (§[3|), the thermal 
energy of a gas particle, ~ kT, where k is Boltzmann's 
constant, is much smaller than the particle's gravita- 
tional energy due to the P1/P2 disk, ~ GMu^l Rd- Here 
G is the gravitational constant, Mo ~ 2 x 10 M a and 
Rd ~ Ipc are the mass and characteristic size of the 
disk, respectively, and fj, is the mean molecular weight of 
the gas. Therefore, as Paczynski (1977) originally envi- 
sioned, gas orbits can be computed in the zero pressure 
limit. 

To calculate Rt, we look for simply closed and non- 
crossing orbits in the combined potential of a central 
SMBH and an eccentric, precessing stellar disk. All or- 
bits are computed in the (x, y) plane of the eccentric stel- 
lar disk. We orient PI and P2 such that PI is centered 
at x = —3 pc (which corresponds to the 2 pc separa- 
tion projected on the sky). The SMBH is fixed at the 
origin. We assume that the P1/P2 disk precesses rigidly 
counterclockwise with pattern speed fl p about the center 
of mass, located at (x cm ,0). In the rotating frame, the 
equations of motion are 

x = - ( ^- + fi^ (x - x cm _) + 2£l p y (1) 



y = -—+n 2 p y-2n p x, (2) 
where the potential $ is given by 

GMrh 

r 

where Mbh is the mass of the central SMBH, r = 
\J ' x 1 + y 2 , and <1>d is the potential due to the eccentric 
disk: 

*D(*>tf) = -G I dx'dy' ^ X ']f]_ (4) 
J \r — r | + h 

where \r — r'\ — wjx — x') 2 + (y — y') 2 , S is the stellar 
mass surface density, and h is the softening length to 
account for the finite thickness of the disk. For now, we 
take a fixed softening length of h = 0.1 pc to cleanly 
demonstrate how the P1/P2 disk can tidally truncates 
a gaseous disk. The P1/P2 disk may be substantially 
thicker, however. At the end of this section, we discuss 
various choices for h and their effects on our results. 

For S, we use the fit to the light distribution from 
PT03 and a mass-to-light ratio of 5.7 (T95; PT03). We 
focus on PT03's non-aligned model as it fits the light 
distribution and kinematic data better for PI and P2. 
In addition, we ignore the contribution to the potential 
from the local bulge because its mass is only 10% of that 
of the disk (< 10 6 M Q ; PT03). We refer the interested 
reader to Appendix [Al for more details. 

The strength of the non-axisymmetric component rel- 
ative to the axisymmetic component depends on the 
mass of the P1/P2 disk relative to that of the SMBH, 
Md/Mbh- PT03 give a stellar mass of w 2 x 10 7 M Q 
for a mass-to-light ratio appropriate for the bulge and 
fit a SMBH mass of lO 8 Af . Salow & Statler (2001, 
2004) construct a series of self-consistent dynamical mod- 
els, which give stellar masses of 1 — 2 x 10 7 M Q and a 
SMBH mass of w 5 x 1O 7 M . B05 give a SMBH mass 
of 1.4 x 10 8 M Q (1.1 - 2.1 x 1O 8 M at the la error level) 
based on the dynamics of P3. The range of mass ratios, 
Md/Mbh> from these different authors is » 0.1 — 0.3. 
We adopt a fiducial ratio of Md/Mbh = 0.1 motivated 
by the stellar disk mass from PT03 and the SMBH mass 
from B05; in this case x cm = —0.07 pc. 

We compute test particle (gas) orbits using equations 
|T]) and (2). The orbit starts at x = —R\ with a ve- 
locity purely in the —y direction (x = 0); see Figure [TJ 
Stars rotate about the SMBH in the counterclockwise 
direction. We take the gas rotation and the P1/P2 disk 
pattern speed to also be in the counterclockwise direc- 
tions. For each computed orbit, we check to see if the 
ending position and velocity are the same as the begin- 
ning position and velocity. We vary the initial velocity y 
until this condition is met. We define i?2 as the position 
where the orbit crosses the positive x-axis. Repeating 
this calculation over a range of R\, we find a family of 
simply closed orbits for a given pattern speed fl p . 

Figure Q] shows such simply closed orbits for 
Md/Mbh = 0.1, M BH = 1O 8 M , and Q p = 
3kms _1 pc _1 . At such a low pattern speed, gas orbits 
are restricted to lie inside a maximal orbit (thick solid 
line), which we denote the Rt orbit. Gas outside this or- 
bit will cross the Rt orbit (as in the outermost two orbits 
in Fig. [1]) or be self-crossing. 
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Fig. 1.— Gas orbits for U p = 3kms _1 pc _1 , M BH = 10 s M , 
Md/Mbh = 0.1, and a softening length of h = 0.1 pc. The largest 
possible orbit, denoted Rt, is shown with the thick solid line. The 
outermost two orbits have their closest approach inside of the Rt 
orbit (they also cross each other). Gas in orbits exterior to the Rt 
orbit will shock and be driven to it. Note that the Rt orbit is very 
eccentric. 
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Fig. 3. — i?2 (periapse distance) as a function of R± (apoapsc 
distance) for Q p = 2 — 9kms~ 1 pc — 1 in 0.5kms~ 1 pc~ 1 intervals, 
for M D /M BH = 0.1, Mbh = 10 s A/ , and h = 0.1 pc. For small 
Q p , R2 obtains a local maximum, which we mark with a solid 
black square. For Q p > 6kms _1 pc _1 , this local maximum does 
not exist. 
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Fig. 2. — Same as Figure [T1 except f2 p = 30kms _1 pc — 1 . Unlike 
in Figure [T] there is no Rt orbit. Gas finds non-crossing orbits at 
all radii. 

Figure [2] shows orbits for a higher pattern speed f2 p = 
30kms _1 pc _1 . These orbits form non-crossing orbits 
spanning the entire disk. Figures Q] and [2] illustrate that 
the nature of gas orbits qualitatively changes when going 
from low to high pattern speeds. Gas occupies restricted 
orbits at low pattern speeds, while at higher pattern 
speeds, gas can span the entire disk. 

We plot i?2 as a function of R\ for different values of 
Q p in Figure [H For small R2 has a local maximum, 
#2,max (marked with solid black squares in Fig. [3]). We 
denote the R\ for which this occurs as i?i, m ax. Beyond 
Ri « 3 pc, i?2 increases with Ri for all f2 p . 

We define the tidal truncation radius as the angle- 
averaged radius of the maximal non- intersecting orbit: 
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Fig. 4. — Rt (solid line) as a function of Q p for M d /M B h = 0.1, 
Mbh = 10 s Mq, and h = 0.1 pc. We also plot Ri, max (dotted 
line), i?2,max (short dashed line), and iit.ana (long dashed line) as 
defined by the crossing condition (eq.[8]). For Q p > 6kms —1 pc -1 , 
gas orbits can persist at all radii. 



J?t = (27r) 1 / R(6)d0. Figure H shows R t as a function 
of f2 p . As ^ increases, Rt increases. The tidal trunca- 
tion radius, Rt, is < 1 pc (similar to the observed max- 
imal radial extent of P3) when fi p < 6kms _1 pc _1 for 
Md/Mbh = 0.1 and M B h = 1O 8 M . For larger pattern 
speeds, Rt does not exist and a gaseous disk can span all 
radii. 

The dependence of Rt on Md, and Mbh can be 
derived using perturbation theory (e.g., Papaloizou & 
Pringle 1977; Binney & Tremaine 1987). In cylindrical 
coordinates (r, cf>), the equations of motion in the frame 
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Fig. 5. — Ratio of $1,0 (rn = 1 component of the potential) to 
$0 (axisymmetric component of the potential) as a function of r 
for' M D /M BH = 0.1. 



rotating at Q p are 



l- 



dr 



„3 



^ ^\^p If 



r 3 r 
20 p rf , 



(5) 
(6) 



where / = r 2 </> is the specific angular momentum. 3 Fol- 
lowing the derivation in Binncy & Trcmaine (1987), 
we take $ -> $ (r) + T,m=i®™{r,<p), r -> r + <5r 
and </> — > + (50 in the epicyclic approximation, where 
$ m = 3> m ,o(0 cos [m(0 — O p )i] are the Fourier compo- 
nents of the potential and O = y / r~ 1 d&o/dr is the or- 
bital frequency at the guiding center, r. If the axisym- 
metric component (3>o) is dominant, the solution for Sr 
is (Binney & Tremaine 1987, their eq. [3.119ab]) 



Sr 



cos (m(0 - fi p )t) 20 
K 2 -m 2 (0-O p ) 2 l^dr + (O - Op); 

$ m ,o(r) + C cos (Kt + tp), (7) 

where k 2 = d 2 §o/dr 2 + 30 2 is the square of the epicyclic 
frequency, C and are constants, and we assume that 
the perturbation is dominated by a single mode m. Note 
that simply closed orbits correspond to C = 0. For the 
eccentric stellar disk, the dominant Fourier mode is m = 
1, whose amplitude we numerically calculate and plot in 
Figure [3 

Orbits first cross their neighbors at pericenter (i.e., 
Sr < 0; see Figure [lj when 



dSr 
dr 



< -1; 



(8) 



i.e., when the epicyclic amplitude grows faster than the 
size of the guiding center orbit (Papaloizou and Pringle 
1977). From our numerically calculated $1.0, we evaluate 
dSr/dr as a function of r and determine an analytic tidal 
truncation radius, i?t.ana, where dSr/dr first equals — 1. 

3 We assume in the following that the disk precession axis passes 
through the origin (SMBH). 



We plot i?t,ana in Figure[4]to compare to the numerically 
calculated Rt- The agreement is good. 

Since the potential is nearly Keplerian (k m O) and 
since Q p /Q < 1, we expand equation to first order 
in fip/f2 and $i/$o- For d/dr ~ 1/r, the condition for 
first orbit crossing (eq.[8]) becomes 



\Sr\ n $1 

r Op $ 



(9) 



where we drop numerical factors. Equation indicates 
that for the fixed R t orbit, O p oc O$i/$ - Since $i/$o 

1 /2 

scales as Md/Mbh and O scales as M B jj , we find: 



0„ oc M 



-1/2 

BH 



Mr 



(10) 



We demonstrated earlier (see Figure 0]) that gas orbits 
are limited to an inner disk similar to the maximal radial 
extent of P3 (i.e., Rt < 1 pc) if the pattern speed O p < 
ekms-ipc" 1 for M BH = 1O 8 M and M D /Af BH = 0.1. 
Inserting the scalings from (fTOj) . we rescale this pattern 
speed to be 

-1/2 

(11) 



o <r ( Mbh 



V 10 7 M G 



km s 1 pc 1 . 



We have confirmed these scalings numerically. For M31, 
if M D w 2 x 10 7 Af o (PT03) and M BH « 1.4 x 10 8 A/ P 



© 



(B05), then equation (fTTj) gives O p < 10kms _1 pc _1 . 

We now return to the issue of the thickness of the 
P1/P2 disk and its impact on R t . The non-aligned model 
of PT03, which we use for the surface density profile, uses 
h/r w 0.4. On the other hand, the disk models of Sa- 
low and Statler (2001, 2004) are razor thin. Bacon et 
al. (2001) also advocate a cold thin disk (h/r ~ 0.1) 
to model PI and P2. Toomre stability arguments give 
a minimum h/r of sa 0.1 (PT03; see also Bacon et al. 
2001). If the P1/P2 disk has persists for 10 10 years, 
two-body relaxation gives a minimum h/r m 0.2 (PT03; 
T95). 

How does varying the softening length affect our re- 
sults? In Figure we plot R t as a function of O p for 
various softening parameters of h/r = 0.1 — 0.4, taking 
Mr> = 2 x 10 7 M Q and M BH = 1-4 x 10 s M q as is appro- 
priate for M31. The results for h = 0.1 pc, which we have 
focused on thus far, differ from those for h/r = 0.1 by 
w 10%. However, as Figure [6] shows, thicker disks show 
more substantial differences. As we increase the soften- 
ing parameter from h/r — 0.1 to 0.4, the maximum O p 
for which R t exists decreases, down to ~ 2kms~ 1 pc~ 1 . 
In addition, the maximum R t also decreases, down to 
k 0.4 pc. To produce R% < 1 pc, similar to the observed 
maximal radius (~ 1 pc) of the P3 disk (B05), we require 
% 3 - 10kms _1 pc for h/r = 0.1 - 0.3. 

Finally, we estimate the time, ffl ow , f° r S as 1° How to Rt 
from larger radii. Compared to the axisymmetric com- 
ponent of the potential, the m = 1 component is smaller 
by of order 10~ 3 — 10~ 2 for r > 1 pc (see Fig. [5|). We 
expect the epicyclic velocity to be Sv ~ d $i/ ( E > o' y orb ~ 
0.1w or b ~ 70kms~ 1 at 1 pc, where w or b is the orbital ve- 
locity. This is much greater than the gas sound speed, 
c s ~ 0.3 — 1 km s _1 at T ~ 10 — 100 K. Gas streams cross 
at supersonic velocities, shock, and dissipate ^i/^o of 
their orbital energy per dynamical time, idyn = O -1 . 
Thus the inflow time ta ow ~ tdyn(^i/^o) \ which is 
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Fig. 6. — -Rt as a function of Q p for parameters appropriate to 
M31 (M D = 2 X 10 7 M Q and A/ BH = 1.4 X 1O 8 M ) for various soft- 
ening parameters of h/r = 0.1 (solid line), 0.2 (long-dashed line), 
0.3 (short-dashed line), and 0.4 (dotted line). We also show the 
softening length, h = 0.1 pc, (dashed-dotted line) for comparison. 
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3. STELLAR MASS LOSS AS THE ORIGIN OF THE GAS 
THAT FORMED P3 

In § 2 we argued that gas at r > Rt is forced down to 
r rj Rt on short timescales, ifl ow ~ 10 6 yrs. Supplied with 
sufficient mass, the gas disk at ~ Rt eventually becomes 
gravitationally unstable, forming stars and giving rise to 
P3. We now explore these questions of gas supply and 
gravitational instability. 

There are a number of potential sources of gas. For 
example, molecular clouds can be gravitationally scat- 
tered onto low angular momentum orbits. The rate of 
such gas supply is difficult to quantify, however. Here we 
focus instead on mass supply due to winds from stars in 
the P1/P2 disk. This source of mass is unavoidable, ex- 
isting even in the absence of external sources. We show 
below that mass loss from the P1/P2 disk is sufficient to 
trigger a starburst having the right properties to explain 
P3. 

The P1/P2 disk consists of w 2 x 10 7 M© of stars with 
an age of - 10 10 yrs (B05, PT03). We compute the stel- 
lar mass loss rate from Starburst99 (Leitherer et al. 1999; 
Vazquez & Leitherer 2005), using the Padova tracks with 
asymptotic giant branch (AGB) stars. A ~ 10 10 yr old 
star cluster of mass 2 x 10 7 M Q loses mass via stellar 
winds at a rate M, w 3x 10~ 5 - 3 x Kr 4 Af© yr -1 . 
The mass loss is primarily due to winds from red gi- 
ants and thermally pulsating AGB stars. There are un- 
certainties in these mass loss rates due to uncertain- 
ties in the properties of the thermally pulsating AGB 
stars. The stellar winds which dominate the mass loss 
have velocities (5 — 25kms _1 for AGB stars; Lamers 
& Cassinelli 1999) much lower than the orbital veloc- 
ity w rb ~ 700 km s -1 . Hence the winds are bound to 
the nuclear region. The winds have typical relative ve- 
locities of ~ (h/r)v OT b ~ 200 km s^ 1 , where h ~ 0.3 r is 
the thickness of the P1/P2 disk. The stellar winds thus 



shock, reaching temperatures of ~ 10 6 [(h/r)/0.3] 2 K. 
The fate of the shocked stellar wind material depends on 
the competition between heating and cooling. For rapid 
cooling, the gas will collapse to form a thin disk. For slow 
cooling, it will heat up because of viscous stresses and 
form a geometrically thick, radiatively inefficient accre- 
tion flow. The gas heats up on a characteristic timescale 
ihcat ~ a^ 1 ^ 1 , where a is the dimensionless viscosity. 4 
The cooling time is given by t coo \ w 3kT/2nA(T) where 
A(T) is the cooling function for an optically thin thermal 
plasma. The density of stellar wind material that accu- 
mulates on a timescale iheat is n w M*iheat / (27rr 3 [/i/r]/i) . 
If ^cooi ^$ ^hoat for gas at this density, it will cool and col- 
lapse to form a thin disk. This condition can be written 
as a constraint on the stellar wind mass loss rate 



M» > M, 



*,crit 



3-irra 2 [h / r] 3 At 2 w 4 rb 
5A(T) ' 



(12) 



Taking a = 0.1, A/bh = 10 8 M Q and r = 2 pc as fiducial 
numbers for M31, we find that M*, cr it ~ 4 x 10~ 7 - 
6 x lO^Moyr" 1 for h/r = 0.2 - 0.4. Since M* w 3 x 
10~ 5 - 3 x 10" 4 M Q yr" 1 for the stars in the P1/P2 disk, 
we conclude that stellar winds from P1/P2 will likely 
cool and collect in a geometrically thin disk. 5 

Cooled gas accumulates at Rt ~ 1 pc until it cither 
accretes onto the SMBH or becomes gravitationally un- 
stable and fragments into stars (see related arguments of 
Nayakshin 2006 and Levin 2007) . For a disk to fragment 
into stars, numerical studies suggest that two conditions 
must be met (Gammie 2001; see also Shlosman & Begel- 
man 1989): 



Q<i, 

^cool 3^dyn? 



(13) 
(14) 



where Q = c s k/7tGE 9 is the Toomre parameter, c s = 
y/kT/fi is the gas sound speed, E ff is the gas surface 
density of the disk, and tdyn = & is the local dynamical 
time. The radiative cooling time, t coo \, is given by 



M^cool 



'IR 



TIR < 1 
TIR > 1 



(15) 



where Tir, = Ki^T, g /2 is the infrared (IR) vertical optical 
depth of the gas disk, T is the midplane temperature of 
the disk, and ktr is the corresponding opacity. The first 
condition (eq. |13| ) is that gas self-gravity must over- 
come rotation and gas pressure. The second condition 
(eq. [T3]) is that cooling is sufficiently rapid to over- 
come heating due to gravitationally induced turbulence 
( "gravitoturbulence" ; e.g., Shlosman et al. 1989; Gam- 
mie 2001). If equation (fT3)) is satisfied, but equation ([TJ| 
is not, then the disk enters a gravitoturbulent state and 
accretes, but does not fragment or form stars. 

4 The heating from the gravitational torques exerted by the ec- 
centric stellar disk is small compared to that from a-viscosity for 
an initially thick disk. 

5 Our critical M* lC rit ' s a factor of > 20 smaller than that usually 
estimated for the transition from a thin to thick disk (e.g., Fig. 3 of 
Menou, Narayan, & Lasota 1999). The latter calculations assume 
h ss r, i.e., that a thick disk is already established, while in our 
problem the stellar winds are initially confined to a region with 
h Ri 0.3 r. Smaller h/r increases n and decreases T at fixed M», 
thus significantly decreasing the cooling time of the gas. 
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The ability of the gas disk to fragment into stars thus 
depends on the heating and cooling of the gas. We con- 
sider two possibilities for the heating: external heating 
by stellar irradiation from the P1/P2 disk and intrinsic 
heating by gravitoturbulence. We take the gas to cool 
radiatively by emission from dust grains. 6 When the gas 
is externally heated by starlight, Q > 1 initially for a suf- 
ficiently low mass disk. The disk mass grows from stellar 
winds until Q ~ 1, when it becomes gravitationally un- 
stable. If equation (TTJJ is also satisfied, then the disk 
fragments. 

When external sources of heat are negligible, gravito- 
turbulence tends to maintain the disk in a marginally 
stable state with Q ~ 1 (Gammie 2001). Initially, the 
disk does not fragment because the cooling time is long 
under these conditions. As the mass in the disk increases 
from stellar winds, the cooling time decreases relative to 
the orbital period, and eventually fragmentation occurs. 
Whether the gas is heated by starlight or by gravitotur- 
bulence, fragmentation is a function of M gas (see also 
Nayakshin 2006 and Levin 2007). 

We first consider stellar irradiation, which is dom- 
inated by the P1/P2 disk. The stars in the P1/P2 
disk are on average Rd ~ 1 — 3 pc from any patch 
of the gaseous disk, whose vertical thickness is <C h, 
the thickness of the stellar disk. For the purposes of 
this section, we will adopt fiducial values of Fd ~ 2 
pc and h/Ru ~ 0.3 motivated by our previous dis- 
cussion in § 2. The flux of starlight incident on 
the disk is F* - (L*/4TrR 2 D )(h/R D ) « 10(L*/3 x 
10 6 L Q )(i? D /2pc)~ 2 ([V J Rr>]/0.3)ergscm- 2 s- 1 , where 
L* is the total stellar luminosity of the P1/P2 disk. 
The disk is easily optically thick to incident starlight for 
typical optical opacities (for dust-to-gas ratios of 0.01; 
Drainc 2003). We define the effective temperature from 
starlight heating as 



F* 



L * 



4vri?2 R D ' 



(16) 



which gives 



20 



3 x 10 6 L f . 



1/A (R D y 1/2 (h/R D v'\ 
teJ [-or) 

The emitted flux is Fir w ctbT min(riR, 1). Equating 
F* with Fir, we find 

r = T ei *min(r m ,l)- 1 / 4 K. (18) 

Note that in the optically thick case, external irradiation 
generates a midplane temperature, T, that is indepen- 
dent of optical depth. At such low temperatures (tens of 
K), the main source of opacity is dust, for which 

2 



Kir = 5.5 



T 



166 K 



cm 2 g _1 T<166K (19) 



(Bell & Lin 1994; Thompson, Quataert, & Murray 2005). 
The disk builds up sufficient mass to fragment when Q ~ 
1, which implies 



M, 



gas, crit 



Mi 



BH 



"orb 



(20) 



6 Gas and dust are thermally well coupled by gas-dust collisions. 
A O.lfim grain equilibrates with surrounding gas in ~ 10 6 s at 
the gas densities n ~ 8 X 10 8 cm -3 and temperatures T ~ 30 K 
characterizing Q ~ 1 disks. 



At this time, the disk may be optically thin or thick. 

We first consider the optically thin case. Combin- 
ing equations (JTHJ), (O, and P0|) , and using S g = 
M gaSiCr it/7rF 2 for a gas disk with radius R, we find 



M gaSjCrit ~5xl0 4 



Mi 



(lOM/. 



6/13 



F 



4/13 



8/13 



20 K 



M . (21) 



From equation (|18p . the corresponding temperature is 



T - 30 



( M BH A" 1713 



V 10 8 M, 



o 



R \ 3/13 f T, 



1 pc 



20 K 



8/13 



K. 



and the cooling time is 

5/13 

^cool ~ 6 



(22) 



V 10 8 M G 



15/13 

ipey Uok 



( m bh Y' 10 ( R V 1J,1 7^, 



-40/13 



yrs. 



The cooling condition (eq. 
Fcooithin where 



(23) 

is satisfied for R > 



R 



cool, thin 



o.i 



/ M 



BH 



V10 8 M Q 



23/69 



F 



20 K 



-80/69 



pc. (24) 



Once the critical gas mass is reached (eq.[2l]) (for R > 
Fcooi,thin), the disk fragments and forms stars. 

For larger gas masses, the disk becomes optically thick. 
Using equations (fl"9]) and ([20]) , we find that the gas mass 
where the optically thin to thick transition occurs is 



M gas , r=1 w 6xl0 4 



R 



4/5 



/ M] 



BH 



-2/5 



M . (25) 



lpcj \W 8 Mq, 

For an optically thick disk, the critical mass for fragmen- 
tation is 

1/2 / R \ 1/2 



M, 



gas, crit 



3 x 10 



4 ( M } 



BH 



V10 8 M G 



0.3 pc 



F, 



20 K 



1/2 



M , (26) 



where the corresponding temperature T = T c ». Note we 
rescaled T c ,* so that the disk is self-consistently optically 
thick. The cooling time is 

~ 3 / M BH 



tcool w 100 



R 



0.3 pc 



V 10 s M 



yrs. 



(27) 



For optically thick cooling in the regime where k oc T 2 , 
the cooling time is independent of the gas mass. In this 
case, the cooling condition (eq.[14|) is satisfied for R > 
F cool , thick = O.2(M BH /lO 8 M )- 1 /3 pc . 

If a disk reaches Q ~ 1, but the cooling condition 
(eq.pT]) is not satisfied, then it cannot fragment im- 
mediately. Instead, gravitoturbulence heats the disk to 
maintain a temperature of 



T « 30 



R 

1 pc 



M 



gas 



/ M 



BH 



5xlO 4 M y \W s Mq 



K 

so as to keep Q ~ 1. To fragment, the disk must ac- 
cumulate additional mass until £lt coo i < 3. Again this 
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cooling may proceed via optically thin or thick emission. 
The critical mass for fragmentation for an optically thin 
disk is 



M, 



gas, cnt 



3 x 10 



4 / M B H 



V 10 8 M G 



11/20 



R 

1 pc 



7/20 



Mr, 



(29) 

For R < i? C ooi, thick, optically thick cooling in the /tocT 2 
regime is too slow for the disk to fragment. Instead, the 
gas mass there will build up and gravitoturbulence will 
heat the disk until k is not longer proportional to T 2 , i.e., 

T > 166 K. Above this temperature, k « 5.5 cm 2 g 1 
is roughly independent of temperature (c.f. Thompson 
et al. 2005), though it varies by factors of a few be- 
tween T ~ 100—1000 K, where the upper bound is set by 
the dust sublimation temperature. Assuming a constant 
opacity, we find the critical gas mass for fragmentation 
to be 

7/8 f R N -5/8 



M. 



gas, cnt 



5 x 10 



4 ( M } 



BH 



V 10 8 Af G 



0.2 pc 



M 0! 
(30) 

for R < i? CO ol, thick- 

We summarize the above results in Figure [7] We use 
the opacity table compiled by Semenov et al. (2003) and 
numerically compute the self-consistent fragmentation 
mass (eqns.[T3] and [14]) with radiative cooling (eq.|15|). 
allowing for starlight (eq. [TH] ) and gravitoturbulent heat- 
ing (eq.[28]l. whichever is greater. Figure [7] shows the 
critical disk mass for gravitational collapse as a function 
of R at T C! » = 5, 20, and 50 K. We also show the i? 8 / 13 
scaling from equation (f2"Tj) , the R}l 2 scaling from equa- 
tion (|26|) . and the i? 7//20 scaling from equation (j29|) in 
their respective regimes. We find that at R w 1 pc, ex- 
ternal irradiation dominates. For R < R CO o\, thick ~ 0.2 
pc, gravitoturbulence heats the central disk tempera- 
ture above which the opacity law, k k T 2 , no longer 
holds. However, we do not recover the scaling suggested 
by equation (|30|) as the opacity is not a constant above 
166 K, but rather varies by factors of a few. 

We have shown that depending on whether disks are 
externally heated (eq. [H]) or internally heated (eq. 
[29]). the fragmentation mass is 3 — 5 x 1O 4 M . Gas 
from stellar mass loss at r > Rt is driven to r ss R t on a 
timescale £ now ~ 10 6 yrs (§2). For mass supply rates of 
M* ~ 10 _4 A/q yr _1 , the steady-state disk mass outside 
Rt is expected to be ~ 100Af Q , well below that required 
to fragment. Thus all of the gas is driven to R w R t ~ 1 
pc, where it collects in a ring. The timescale for gas to 
viscously spread once it accumulates at Rt is 

2 _a / Af BH 

i 



idyn 

a 



— « 5000a" 



Rt 
1 pc 



1/2 



V io 8 m q 
i 



1/2 



T 



30 K 



Myrs, (31) 



where a < 1 is the dimensionless viscosity. We compare 
this to the time needed to accumulate a gravitationally 



unstable disk, t a 



-Mgas,crit/-M*- The mass loss rate 



varies as the stellar population ages. Using Starburst99, 
we find that as a 2 x 10 7 M Q stellar cluster ages from 
3 x 10 9 to 10 10 yrs, the mass loss rate, M*, ranges from 
3 x 10~ 4 to 3 x lQ- 5 M Q yr _1 . The range in M* and the 



10 5 





/. 


- 








- 








A/ 






1 


i 


1 




R [pc] 



10 4 



Fig. 7. — Critical gas mass for fragmentation as a function of 
disk radius R at T Cl » = 5 (lower dashed line), 20 (thick solid line), 
and 50 K (upper dashed line) for M BH = 10 s A/ Q . At _R < 0.2 pc, 
local accretion heating dominates, while at larger radii, irradiation 
dominates. The dotted line shows the disk mass at which "tir = 
1. We also show the _R 8 / 13 scaling from equation H21I I. the R 1 / 2 
scaling from equation 1261 1. and the _R 7 ^ 20 scaling from equation 



(29). 

range in M gas , cri t yield t 

a 

the ratio 



100 - 2000 Myrs. Hence 



0.1a 



Af. 



( T 



10- 4 M e yr- 



V30K 



3/2 



(32) 



where we have used the critical mass for Q ~ 1 
from equation (|20p . Thus, even for a ~ 1, the ratio 
^accumAvisc % 0.3, which implies that gas fragments at 
~ Rt before accreting. 

The current mass in stars in P3 is estimated to be 
w 4200M Q (B05). The disk mass required for gravi- 
tational instability is ~ 6 — 10 times higher than this 
(Fig. [7]). Such a difference suggests a star formation ef- 
ficiency ~ 10 — 20%. A larger star formation efficiency 
can be accommodated if the initial mass function (IMF) 
of stars in this extreme environment is top heavy, which 
would help reconcile our estimated M gaSjCr it with the cur- 
rent inferred mass of P3. There is theoretical (Nayakshin 
2006; Levin 2007) and observational (Nayakshin & Sun- 
yaev 2005; Paumard et al. 2006) evidence for a top heavy 
IMF in the GC starburst. 

Because fd yn <C iaccum, it is also possible that only a 
small fraction of the gas fragments into stars once the 
disk mass exceeds Af gaSi crit, leaving behind a stable disk 
with a mass only modestly below M gas , crit- The excess 
gas would remain in the nuclear region, given the severe 
difficulty that stellar winds and supernovae would have 
in removing it from so deep in the gravitational potential 
of the black hole. In this case, the gas disk in the nucleus 
of M31 would maintain a mass of ~ M gaS)Cr jt. 

Since the starburst that produced P3 occurred, we ex- 
pect M* x 200 Myr ~ 10 4 Af Q of gas to have accumulated 
near R t ~ 1 pc from ongoing stellar mass loss. This 
molecular gas would be analogous to the circumnuclear 
disk (CND) in the GC, with temperatures of T « 30 K 
and extremely high densities of n ~ 10 9 cm -3 (Q > 1). 
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This gas would be bright in CO and HCN; the CO (1-0) 
flux would be w 2 mJy (for an optically thick line). 

4. DISCUSSION AND CONCLUSIONS 

We have argued that the origin of the young stars 
in P3 in M31 is rooted in the P1/P2 disk. The 
non-axisymmetric component of the potential from the 
P1/P2 disk restricts non- intersecting gas orbits to dis- 
tances r < i? t ~ 1 pc from the central SMBH if the 
pattern speed (i.e., precession frequency) of the P1/P2 
disk is fip < 3 — 10kms _1 pc^ 1 . At larger radii, gas finds 
itself in intersecting orbits; it shocks, loses energy, and 
gets driven to Rt ~ 1 pc. This is comparable to the 
maximum radial extent of the A stars of P3. 

Stellar mass loss from the P1/P2 disk can supply the 
gas that formed P3. Stellar winds supply mass at a rate 
of ~ lO _4 Af yr" 1 for a ~ 10 10 yr old population. This 
gas accumulates in a disk at r < Rt- The conditions 
for fragmentation (eqns.[l3] and [Mj) give a critical gas 
mass of ~ 5 x 1O 4 M . Hence, every ~ 500 Myr, the 
disk accumulates enough mass to fragment and produces 
a starburst. This recurrence time is consistent with the 
age of the A stars of 200 Myr (B05). In addition, the 
observed alignment of the P3 disk with the P1/P2 disk 
is consistent with our argument that the P1/P2 disk sup- 
plies the gas which formed P3. 

Several predictions arise naturally out of this model. 
First, the pattern speed of the eccentric stellar disk 
should be H p < 3 — 10kms _1 pc _1 . Current observa- 
tional constraints on the pattern speed are weak (see 
Appendix A). Second, there should be ~ 1O 4 M of gas 
in the nucleus of M31 from accumulation of gas since 
the last starburst that produced P3. This molecular gas 
would be analogous to the CND at the GC. Such gas 
would have a temperature of ~ 30 K, be at extremely 
high densities n ~ 10 9 cm -3 , and be bright in CO and 
HCN with a CO (1-0) flux of w 2 mJy and a line width 
of w 1000 km s~ . Under these conditions, dust will emit 
with a flux of ~ 20 mJy at 70 /im, ~ 100 mJy at 160 /im, 
and ~ 0.1 mJy at 1 mm. Finally, starbursts at < 1 pc in 
M31 should occur every 10 s — 10 9 yrs. These older gen- 
erations of stars and stellar remnants may be detectable. 

Another interesting possibility, which was suggested 
by the referee, is that the P3 stars may be chemically 
anomalous because they are constructed from the re- 
cycled AGB winds of the P1/P2 disk. A similar self- 
pollution scenario has been proposed in the context of 
globular clusters (Cottrell & Da Costa 1981; Ventura et 
al. 2001). The composition of the AGB wind depends 
strongly on the initial mass of the star and the amount 
of dredge-up between the core and the envelope (Fenner 
et al. 2004), which is uncertain. In light of these uncer- 
tainties, it is interesting to note that HST observations 
of omega Centauri, the largest Galactic globular cluster, 
show a double main sequence (Bedin et al. 2004), sug- 
gesting two episodes of star formation. The blue main 
sequence, which arises from the second episode of star 
formation, may possess a considerable enhancement of 
helium (Bekki & Norris 2006; Karakas et al. 2006), 
which could come from the AGB winds of the red main 
sequence, i.e., the stars of the first episode of star forma- 
tion. It would be interesting to search for an analogous 
chemical anomaly in the P3 stars of M31. 

Observations indicate that P3 appears to be a circular 



disk around the SMBH. Gas orbits around the SMBH 
are eccentric in the presence of the P1/P2 disk as illus- 
trated in Figure [1] However, stars that form from an 
eccentric gas disk may not have an eccentric distribution 
themselves at later times. Once the Q ~ 1 gas disk turns 
into stars, these stars will precess differentially because 
of their own self-gravity. We estimate the differential 
precession rate to be 



dR \M BH J \ R J 



(33) 



where w is the precession frequency of a star at radius 
R, and AR is the initial radial extent of the P3 stel- 
lar disk with mass Mp 3 . A spread of AR/R ~ 0.1 in the 
orbits of the P3 stars can be generated by viscous spread- 
ing of the deposited gas prior to star formation. Taking 
Mp3/Mbh ~ 10~ 4 , we find that stars in P3 differentially 
precess out of their initially apsidally aligned structure 
over a timescale 10 5 f2 -1 ~ 10 8 years, comparable to the 
age of the A stars. 

Over 10 Gyrs, a mass loss rate of 3 x 10~ 5 to 3 x 
lO _4 M yr^ 1 from the P1/P2 disk will redistribute w 
3 x 10 5 - 3 x 1O 6 M of mass to the P3 disk, which is of 
order 10% the mass of the P1/P2 disk. If a large fraction 
of this mass is retained, it may affect the eccentricity of 
the P1/P2 disk. The backreaction of a growing P3 disk 
on the eccentric P1/P2 disk is beyond the scope of this 
paper, but is an interesting question for further study. 

Our model may be applicable to other galaxies with 
double nuclei in addition to M31. Lauer et al. (1996) 
observed a double nucleus in NGC 4486B with a 12 pc 
separation. Debattista et al. (2006) detected a double 
nucleus in VCC 128 with a 32 pc separation. Thatte et 
al. (2000) also detected a double nucleus in the star- 
burst galaxy M83 with a 5.4 pc separation. If these dou- 
ble nuclei are associated with non-axisymmetric stellar 
distributions, very compact nuclear starbursts and dense 
nuclear molecular gas may be common features of galac- 
tic nuclei. 

Finally, we briefly discuss our model in the context of 
the GC. Observations suggest that the 1.3 x 1O 4 M of 
young massive stars in the GC are concentrated between 
r ~ 0.04 - 0.4 pc (Ghez et al. 2005; Paumard et al. 
2006), similar in mass and radial extent to M31. A non- 
axisymmetric component of the potential may explain 
the radial extent of these young stars, which otherwise 
can only be accounted for by the assumption that gas 
is supplied on very low angular momentum orbits. If 
the non-axisymmetric component were due to an eccen- 
tric disk of old stars, it would likely remain undetected 
because of extinction. 

We thank R. Genzel, A. Loeb, B. Paczynski, L. 
Strubbe, and S. Tremaine for useful discussions. We 
thank L. Blitz and A. Bolatto for useful discussions and 
for performing CO observations on M31. We thank the 
anonymous referee for useful comments. We thank B. 
Johnson for presenting a talk on M31 in the Astro Read- 
ing Group (ARG) that led to this project. We also thank 
G. Bower, G. Howes, B. Metzger, and B. Schmekel for 
leading other seminars in the ARG. We would also like to 
acknowledge all ARG participants. P.C. thanks the In- 
stitute for Advanced Study and the Canadian Institute 



Origin of Young Stars in M31 



9 



TABLE 1 

Model parameters from Peiris h Tremaine (2003). 



Parameter 


Aligned Model 


Unaligned Model 


a 


0.288 


0.197 


a e [pc] 


3.97 


4.45 


a g [pc] 


1.51 


1.71 


w [pc] 


1.53 


1.52 


a [pc] 


4.61 


1.37 


C2 [pel 


3.79 


4.24 
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APPENDIX 
DYNAMICS OF PI AND P2 

The most widely accepted model for the double nucleus of M31 is that of an eccentric stellar disk (T95). T95 fits the 
light distribution with three apsidally aligned ellipses with the SMBH at one focus. Stars on these elliptical Keplerian 
orbits pile up at apoapse, giving rise to PI. A suitable nesting of orbits gives rise to P2 at periapse. PT03 performed 
a more careful calculation to fit the light and velocity distributions, using Keplerian ellipses again. These models do 
not take disk self-gravity into account. Since the stellar disk mass is nearly 2 x 10 7 M©, comparable to the SMBH mass 
of ~ lO 8 Af0 (B05), disk self-gravity has significant effects on the dynamics of both stars and gas. 

Other workers have developed self-consistent models for the dynamics of such an eccentric stellar disk. Statler (1999) 
created self-consistent models for a stellar disk assuming that stars follow periodic orbits. This work was expanded by 
Salow & Statler (2004), who find high pattern speeds f2 p ~ 30kms~ 1 pc~ 1 . Sambhus & Sridhar (2002) used families 

of prograde and retrograde loop orbits to fit the light distribution and found Q p « 16kms~ . Bacon et al. (2001) 
showed through N-body simulations that the m = 1 perturbation can be long-lived, ~ 100 Myrs for low pattern speeds 
fi p < 3kms _1 pc _1 

pattern speed of Q p /Q,q 



. Similarly, Jacobs & Sellwood (2001) found long-lived m — 1 modes in an eccentric ring with a 
0AM D /M B a, where fl 



of ftp as 16kms 



\J GM/Rq and Rq ~ 2 pc. This corresponds to a pattern speed 
pc^ 1 for the parameters appropriate to M31. Finally, Sambhus & Sridhar (2000) utilized a variant 
of the Tremaine- Weinberg method (Tremaine & Weinberg 1984) to obtain an observational constraint on the pattern 
speed of < 30kms _1 pc _1 (for a disk inclination relative to the sky plane of 77 degrees). 

As the constraints on the pattern speed are weak, we have elected to take PT03's fits to the P1/P2 disk and to 
experiment with a range of values of Q p . We use their formulae for the disk eccentricity profile and surface density 
(their eqs.[12] and [17]): 



j (a) = a (a e — a) exp 



2w 2 



S(a) = S 



a 2 exp (— a/ao) 
l + exp[4(a-c 2 )]' 



(Al) 



(A2) 



where a is the semi- major axis in pc, e m is the mean eccentricity, and the effective surface density, S(a), is defined as 
dM = 2iraY>(a)da, where dM is the mass between a and a + da. We take the maximum semi-major axis to be a = 8 pc. 
The normalization So is set by Md/Mbh> which we vary. For convenience we list the fitting parameters (their Table 
2) in our Table 1. We focus on PT03's nonaligned model in this paper, but we have studied their aligned model as 
well. In the aligned model, for R t < 1 pc, we find the required pattern speed is il p < 6kms _1 pc~ 1 for Md/Msh = 0.1 
and Mbh = 10 8 M Q for the softening length of h = 0.1 pc, similar to the nonaligned model. 
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